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Abstract: We propose a semiclassical method to calculate 5-matrix elements for two- 
stage gravitational transitions involving matter collapse into a black hole and evaporation 
of the latter. The method consistently incorporates back-reaction of the collapsing and 
emitted quanta on the metric. We illustrate the method in several toy models describing 
spherical self-gravitating shells in asymptotically flat and AdS space-times. We find that 
electrically neutral shells reflect via the above collapse-evaporation process with probability 
exp(— B), where B is the Bekenstein-Hawking entropy of the intermediate black hole. This 
is consistent with interpretation of exp(B) as the number of black hole states. The same 
expression for the probability is obtained in the case of charged shells if one takes into 
account instability of the Cauchy horizon of the intermediate Reissner-Nordstrom black 
hole. Our semiclassical method opens a new systematic approach to the gravitational 
5-matrix in the non-perturbative regime. 
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1 Introduction 

Gravitational scattering has been a subject of intensive research over several decades, see [1] 
and references therein. Besides being of its own value, this study presents an important 
step towards resolution of the information paradox — an apparent clash between unitarity 
of quantum evolution and black hole thermodynamics [2, 3]. Recently the interest in this 
problem has been spurred by the AMPS (or “firewall”) argument [4, 5] which suggests that 
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Figure 1. Complete gravitational transition involving formation and evaporation of a black hole. 


reconciliation of the black hole evaporation with unitarity would require drastic departures 
from the classical geometry in the vicinity of an old black hole horizon (see [6-8] for related 
works). Even the minimal versions of such departures appear to be at odds with the 
equivalence principle. This calls for putting all steps in the logic leading to this result on 
a firmer footing. 

Unitarity of quantum gravity is strongly supported by the arguments based on 
AdS/CFT correspondence [9, 10]. This reasoning is, however, indirect and one would 
like to develop an explicit framework for testing unitarity of black hole evaporation. In 
particular, one would like to see how the self-consistent quantum evolution leads to the 
thermal properties of the Hawking radiation and to test the hypothesis [11] that the infor¬ 
mation about the initial state producing the black hole is imprinted in subtle correlations 
between the Hawking quanta. 

A natural approach is to view the formation of a black hole and its evaporation as 
a two-stage scattering transition, see Fig. 1. The initial and final states \Eq and Tj of 
this process represent free matter particles and free Hawking quanta in flat space-time. 
They are the asymptotic states of quantum gravity related by an 5-matrix [1, 12, 13]. 
Importantly, the black hole itself, being metastable, does not correspond to an asymtotic 
state. The 5-matrix is unitary if black hole formation does not lead to information loss. 
The importance of collapse stage for addressing the information paradox was emphasized 
in [14-18], 

However, calculation of the scattering amplitude for the process in Fig. 1 encounters a 
formidable obstacle: gravitational interaction becomes strong in the regime of interest and 
the standard perturbative methods break down [19, 20]. General considerations [20] sup¬ 
ported by the perturbative calculations [19] show that scattering of two trans-Planckian 
particles is accompanied by an increasingly intensive emission of soft quanta as the thresh¬ 
old of black hole formation is approached. While this is consistent with the qualitative 
properties of the Hawking radiation dominated by many quanta, a detailed comparison is 
not available. A new perturbative scheme adapted to processes with many particles in the 
final state has been recently proposed in [21]; however, its domain of applicability is yet to 
be understood. 

In this paper we follow a different route. We propose to focus on exclusive processes 
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where both initial and final states contain a large number of soft particles. Specifically, one 
can take 'kj and \k/ to be coherent states with large occupation numbers corresponding to 
the semiclassical wavepackets. We assume that the total energy of the process exceeds the 
Planck scale, so that the intermediate black hole has mass well above Planckian. Then the 
overall process is expected to be described within the low-energy gravity. Its amplitude 
can be evaluated using the semiclassical methods and will yield the black hole 5-matrix 
in the coherent-state basis [22]. A priori, we cannot claim to describe semiclassically the 
dominant scattering channel with Hawking-like final state which is characterized by low 
occupation numbers 1 . Still, it seems a safe bet to expect that, within its domain of validity 2 , 
this approach will provide valuable information on the black hole-mediated amplitudes. 

A crucial point in the application of semiclassical methods is the correct choice of the 
semiclassical solutions. Consider the amplitude 

= J V<S>iV<5> f ^[$/] J V<S>e ism/h ( 1 . 1 ) 

of transition between the initial and final asymptotic states with wave functionals dq [<!>,:] 
and T/[<!>/]. The path integral in Eq. (1.1) runs over all fields $ of the theory including 
matter fields, metrics and ghosts from gauge-fixing of the diffeomorphism invariance 3 4 ; S 
is the action. In the asymptotic past and future the configurations d? in Eq. (1.1) must 
describe collections of free particles in flat space-time. In the semiclassical approach one 
evaluates the path integral (1.1) in the saddle-point approximation. The saddle-point 
configuration must inherit the correct flat-space asymptotics and, in addition, extremize S 
i.e. solve the classical equations of motion. However, a naive choice of the solution fails to 
satisfy the former requirement. Take, for instance, the solution <f> c ; describing the classical 
collapse. It starts from collapsing particles in flat space-time, but arrives to a black hole 
in the asymptotic future. It misses the second stage of the scattering process — the black 
hole decay. Thus, it is not admissible as the saddle point of (1.1). One faces the task of 
enforcing the correct asymptotics on the saddle-point configurations. 

Furthermore, the saddle-point solutions describing exclusive transitions are generically 
complex-valued and should be considered in complexified space-time [22, 23] . 1 Even sub¬ 
ject to appropriate boundary conditions, such solutions are typically not unique. Not all 
of them are relevant: some describe subdominant processes, some other, when substituted 
into the semiclassical expression for the amplitude, give nonsensical results implying expo¬ 
nentially large scattering probability. Choosing the dominant physical solution presents a 
non-trivial challenge. 

4 It is conceivable, however, that this dominant amplitude can be obtained with a suitable limiting 
procedure. A hint comes from field theory in flat space where the cross section of the process 2 —> many is 
recovered from the limit of the cross section many —> many [23-25]. 

2 In line with the common practice, the applicability of the semiclassical method will be verified a 
posteriori by subjecting the results to various consistency checks. 

3 The precise definition of the functional measure is not required in the leading semiclassical approxima¬ 
tion, which we will focus on. 

4 This property is characteristic of dynamical tunneling phenomena which have been extensively studied 
in quantum mechanics with multiple degrees of freedom [26]. 
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The method to overcome the two above problems has been developed in Refs. [27-32] 
in the context of scattering in quantum mechanics with multiple degrees of freedom; it 
was applied to field theory in [33, 34], In this paper we adapt this method to the case of 
gravitational scattering. 

It is worth emphasizing the important difference between our approach and perturba¬ 
tive expansion in the classical black hole (or collapsing) geometry, which is often identified 
with the semiclassical approximation in black hole physics. In the latter case the evapo¬ 
ration is accounted for only at the one-loop level. This approach is likely to suffer from 
ambiguities associated with the separation of the system into a classical background and 
quantum fluctuations. Instead, in our method the semiclassical solutions by construction 
encapsulate black hole decay in the leading order of the semiclassical expansion. They 
consistently take into account backreaction of the collapsing and emitted matter quanta 
on the metric. Besides, we will find that the solutions describing the process of Fig. 1 are 
complex-valued. They bypass, via the complexified evolution, the high-curvature region 
near the singularity of the intermediate black hole. Thus, one does not encounter the 
problem of resolving the singularity. 

One the other hand, the complex-valued saddle-point configurations do not admit a 
straightforward interpretation as classical geometries. In particular, they are meaningless 
for an observer falling into the black hole: the latter measures local correlation functions 
given by the path integrals in the in-in formalism — with different boundary conditions and 
different saddle-point configurations as compared to those in Eq. (1.1). This distinction 
lies at the heart of the black hole complementarity principle [35]. 

Our approach can be applied to any gravitational system with no symmetry restric¬ 
tions. However, the task of solving nonlinear saddle-point equations is rather challenging. 
In this paper we illustrate the method in several exactly tractable toy models describing 
spherical gravitating dust shells. We consider neutral and charged shells in asymptotically 
flat and anti-de Sitter (AdS) space-times. Applications to field theory that are of primary 
interest are postponed to future. 

Although the shell models involve only one collective degree of freedom — the shell 
radius — they are believed to capture some important features of quantum gravity [36- 
39]. Indeed, one can crudely regard thin shells as narrow wavepackets of an underlying 
field theory. In Refs. [39-41] emission of Hawking quanta by a black hole is modeled as 
tunneling of spherical shells from under the horizon. The respective emission probability 
includes back-reaction of the shell on geometry, 

( 1 . 2 ) 

where Bi and Bf are the Bekenstein-Hawking entropies of the black hole before and af¬ 
ter the emission. It has been argued in [42] that this formula is consistent with unitary 
evolution. 

In the context of shell models we consider scattering processes similar to those in 
Fig. 1: a classical contracting shell forms a black hole and the latter completely decays due 
to quantum fluctuations into an expanding shell. The initial and final states and ^f of 
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the process describe free shells in flat or AdS space-times. Our result for the semiclassical 
amplitude (1.1) has the form 

(^/|<S|^i) ~ e lSre3 ^ h . (1.3) 

We stress that it includes backreaction effects. The probability of transition is 

T*fi — e-2Im S reg/fo 

We show that for neutral shells it reproduces Eq. (1.2) with Bi equal to the entropy of the 
intermediate black hole and Br = 0. This probability is exponentially small at M > 1 
when the semiclassical approximation is valid. It is consistent with the result of Refs. [36— 
39] since the first stage of the process, i.e. formation of the intermediate black hole, proceeds 
classically. 

For charged black holes the same result is recovered once we take into account instabil¬ 
ity of the inner Cauchy horizon of the intermediate Reissner-Nordstrom black hole [43-49]. 
Our results are therefore consistent with the interpretation of Hawking radiation as tun¬ 
neling. However, we obtain important additional information: the phases of the 5-matrix 
elements which explicitly depend, besides the properties of the intermediate black hole, on 
the initial and final states of the process. 

The paper is organized as follows. In Sec. 2 we introduce general semiclassical method 
to compute 5-matrix elements for scattering via black hole formation and evaporation. In 
Sec. 3 we apply the method to transitions of a neutral shell in asymptotically flat space- 
time. We also discuss relation of the scattering processes to the standard thermal radiation 
of a black hole. This analysis is generalized in Sec. 4 to a neutral shell in asymptotically 
AdS space-time where scattering of the shell admits an AdS/CFT interpretation. A model 
with an electrically charged shell is studied in Sec. 5. Section 6 is devoted to conclusions 
and discussion of future directions. Appendices contain technical details. 

2 Modified semiclassical method 

2.1 Semiclassical 5—matrix for gravitational scattering 

The 5-matrix is defined as 


('k/l^l'Ih) = lim t . ('&f\U 0 (0,tf) U(tf,ti) Uo(ti,0)\'&i) , (2.1) 

t f —y -(- oo 


where U is the evolution operator; free evolution operators Uq on both sides transform 
from Schrodinger to the interaction picture. In our case U describes quantum transition 
in Fig. 1, while JJq generates evolution of free matter particles and Hawking quanta in the 
initial and final states. The time variable t E [U, tf] is chosen to coincide with the time of 
an asymptotic observer at infinity. 

Using the path integrals for the evolution operators and taking their convolutions 
with the wave functionals of the initial and final states, one obtains the path integral 
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Figure 2. The contour used in the calculation of the 5-matrix elements. Quantum transition from 
ti to tf is preceded and followed by the free evolution. 


representation for the amplitude 5 (2.1), 




j ,ti)+iS o (tf,0+) [$ + ] 


f V$e iStotm , 


( 2 . 2 ) 


where <h = {<p, g^u} collectively denotes matter and gravitational fields along the time 
contour in Fig. 2. The quantum measure should include some ultraviolet regularization 
as well as gauge-fixing of the diff-invariance and respective ghosts. A non-perturbative 
definition of this measure presents a well-known challenge. Fortunately, the details of 
D<F are irrelevant for our leading-order semiclassical calculations. The interacting and 
free actions S and So describe evolution along different parts of the contour. The initial- 
and final-state wave functionals T* and Tj depend on the fields <Fip = <b(f = CR) at the 
endpoints of the contour. In the second equality of Eq. (2.2) we combined all factors in the 
integrand into the “total action” Sfot ['!’]• Below we mostly focus on nonlinear evolution 
from ti to tf and take into account contributions from the dashed parts of the contour in 
Fig. 2 at the end of the calculation. 

To distinguish between different scattering regimes, we introduce a parameter P char¬ 
acterizing the initial state [50] — say, its average energy. If P is small, the gravitational 
interaction is weak and the particles scatter trivially without forming a black hole. In 
this regime the integral in Eq. (2.2) is saturated by the saddle-point configuration <F C ; 
satisfying the classical field equations with boundary conditions related to the initial and 
final states [22], However, if P exceeds a certain critical value P*, the classical solution 
<h c ; corresponds to formation of a black hole. It therefore fails to interpolate towards the 
asymptotic out-state \k f living in flat space-time. This marks a breakdown of the standard 
semiclassical method for the amplitude (2.2). 

To deal with this obstacle, we introduce a constraint in the path integral which ex¬ 
plicitly guarantees that all field configurations from the integration domain have flat 
space-time asymptotics [30-32], Namely, we introduce a functional Ti nt [&\ with the fol¬ 
lowing properties: it is (i) diff-invariant; (ii) positive-definite if <F is real; (iii) finite if <f> 
approaches flat space-time at t —> Too; (iv) divergent for any configuration containing a 
black hole in the asymptotic future. Roughly speaking, Tj n t [<h] measures the “lifetime” of 
a black hole in the configuration <h. Possible choices of this functional will be discussed in 


5 From now on we work in the Planck units h = c = Gn = kB = 1- 
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the next subsection; for now let us assume that it exists. Then we consider the identity 

+oo +oo zoo 

1 = J dT 0 S(T mt [<P}-T 0 ) = j (ITq J A e <T 0 - Trnt m ) (2.3) 

0 0 —zoo 

where in the second equality we used the Fourier representation of the 5-function. Inserting 
Eq. (2.3) into the integral (2.2) and changing the order of integration, we obtain, 

= J d Y^ e eT ° J T>4> . (2.4) 

The inner integral over <f> in Eq. (2.4) has the same form as the original path integral, but 
with the modified action 

S e [4>] = Stotm + ieT int [^] . (2.5) 

By construction, this integral is restricted to configurations <f> with Tj n t[<b] = Tq < oo, i.e. 
the ones approaching flat space-time in the asymptotic past and future. In what follows 
we make this property explicit by deforming the contour of e-integration to Re e > 0. Then 
configurations with black holes in the final state i.e. with Tmt[<f>] = +oo, do not contribute 
into the integral at all. 

By now, we have identically rewritten the integral (2.2) in the form (2.4). Its value 
clearly does not depend on the form of the regulating functional T) nt [4>]. 

Now the amplitude (2.4) is computed by evaluating the integrals over 4>, e and To one 
after another in the saddle-point approximation. The saddle-point configuration 4> e of the 
inner integral extremizes the modified action (2.5), while the saddle-point equation for e 
gives the constraint 

Tint[^e] = To . (2.6) 

This implies that has correct flat-space asymptotics. The integral over To is saturated 
at e = 0. Importantly, we do not substitute e = 0 into the saddle-point equations for 4> e , 
since in that case we would recover the original classical equations together with incorrect 
asymptotics of the saddle-point solutions. Instead, we understand this equation as the 
limit 

e -> +0 (2.7) 

that must be taken at the last stage of the calculation. The condition Re e > 0 is required 
for convergence of the path integral (2.4). We obtain the saddle-point expression (1.3) for 
the amplitude with the exponent 1 ’ 


S reg = ^lmig Stot [3>e] , (2.8) 

where the limit is taken in the end of the calculation. 

Our modified semiclassical method addresses another problem mentioned in the Intro¬ 
duction: it allows to select the relevant saddle-point configurations from the discrete set of 

’’Below we consider only the leading semiclassical exponent. The prefactor in the modified semiclassical 
approach was discussed in [30-32]. 


- 7 - 




complex-valued solutions to the semiclassical equations. As discussed in [27-29], at e > 0 
the physical solutions can be obtained from the real classical ones by analytic continuation 
in the parameters of the in- and out-states. This suggests the following strategy. We 
pick up a real classical solution 4>o describing scattering at a small value of the parameter 
P < P*. By construction, <3?o approaches flat space-time at t —> =Foo and gives the domi¬ 
nant contribution to the integral (2.4). Next, we modify the action and gradually increase 
e from e = 0 to the positive values constructing a continuous branch of modified solutions 
<f> e . At e —> +0 these solutions reduce to 4>o and therefore saturate the integral (2.4). 
We increase the value of P to P > P* assuming that continuously deformed saddle-point 
configurations remain physical'. In this way we obtain the modified solutions and the 
semiclassical amplitude at any P. We stress that our continuation procedure cannot be 
performed with the original classical solutions which, if continued to P > P*, describe 
formation of black holes. On the contrary, the modified solutions 4*^ interpolate between 
the flat-space asymptotics at any P. They are notably different from the real classical so¬ 
lutions at P > P*. At the last step one evaluates the action Stot on the modified solutions 
and sends e —> +0 obtaining the leading semiclassical exponent of the 5-matrix element, 
see Eqs. (1.3), (2.8). 

2.2 The functional T) n t[4>] 

Let us construct the appropriate functional P) n t [$]. This is particularly simple in the case 
of reduced models with spherically—symmetric gravitational and matter fields. The general 
spherically-symmetric metric has the form 

ds 2 = g ab{y)dy a dy b + r 2 (y)dQ 2 , (2.9) 

where d£l is the line element on a unit two-sphere and g a b is the metric in the transverse 
two-dimensional space * * * 8 . Importantly, the radius r(y ) of the sphere transforms as a scalar 
under the diffeomorphisms of the y-manifold. Therefore the functional 

Tint = J d 2 y y/^gw(r)F(g ab d a rd b r) , (2.10) 

is diff-invariant. Here w(r ) and F( A) are non-negative functions, so that the functional 
(2.10) is positive-definite. We further require that F( A) vanishes if and only if A = 1. 
Finally, we assume that w(r) significantly differs from zero only at r < r w , where r w is 
some fixed value, and falls off sufficiently fast at large r. An example of functions satisfying 
these conditions is w(r) = 5(r — r w ), F( A) = (A — l) 2 . 

To understand the properties of the functional (2.10), we consider the Schwarzschild 
frame where r is the spatial coordinate and the metric is diagonal. The functional (2.10) 

'In other words, we assume that no Stokes lines [51] are crossed in the course of deformation. This 

conjecture has been verified in multidimensional quantum mechanics by direct comparison of semiclassical 

and exact results [27-32, 52]. 

8 We use the signature (—, +,...) for the metrics g and g a b- The Greek indices g, v, ... are used for the 
four-dimensional tensors, while the Latin ones a, b,... = 0, 1 are reserved for the two-dimensional space of 
the spherically reduced model. 
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takes the form. 


Tint = J dtdr^gw(r)F(g n ) . (2.11) 

Due to fast falloff of w(r) at infinity the integral over r in this expression is finite. However, 
convergence of the time integral depends on the asymptotics of the metrics in the past and 
future. In flat space-time g 11 = 1 and the integrand in Eq. (2.10) vanishes. Thus, the 
integral over t is finite if g a /, approaches the flat metric at t —> ±oo. Otherwise the integral 
diverges. In particular, any classical solution with a black hole in the final state leads 
to linear divergence at t — > +oo because the Schwarzschild metric is static and g 11 ^ 1. 
Roughly speaking, X) n t can be regarded as the Schwarzschild time during which matter 
fields efficiently interact with gravity inside the region r < r w . If matter leaves this region 
in finite time, Tj n j takes finite values. It diverges otherwise. Since the functional (2.10) is 
diff-invariant, these properties do not depend on the particular choice of the coordinate 
system. 

The above construction will be sufficient for the purposes of the present paper. Beyond 
the spherical symmetry one can use the functionals Ti n t[<F\ that involve, e.g., an integral 
of the square of the Riemann tensor, or the Arnowitt-Deser-Misner (ADM) mass inside a 
large volume. Recall the final result for the 5-matrix does not depend on the precise choice 
of the functional T n t. Of course, this functional should satisfy the conditions (i)-(iv) listed 
before Eq. (2.3). 

3 Neutral shell in flat space—time 
3.1 The simplest shell model 

We illustrate the method of Sec. 2 in the spherically symmetric model of gravity with thin 
dust shell for matter. The latter is parameterized by a single collective coordinate — the 
shell radius r(r) — depending on the proper time along the shell r. This is a dramatic 
simplification as compared to the realistic case of matter described by dynamical fields. 
Still, one can interprete the shell as a toy model for the evolution of narrow wavepackets 
in held theory. In particular, one expects that the shell model captures essential features 
of gravitational transition between such wavepackets. 9 

The minimal action for a spherical dust shell is 

Sshell = ~m J dr (3.1) 

where m is the shell mass. However, such a shell always collapses into a black hole and 
hence is not sufficient for our purposes. Indeed, as explained in Sec. 2.1, in order to select 
the physically relevant semiclassical solutions we need a parameter P such that an initially 
contracting shell reflects classically at P < P* and forms a black hole at P > P*. We 
therefore generalize the model (3.1). To this end we assume that the shell is assembled 
from particles with nonzero angular momenta. At each point on the shell the velocities of 

9 Note that our approach does not require complete solution of the quantum shell model which may be 
ambiguous. Rather, we look for complex solutions of the classical equations saturating the path integral. 
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the constituent particles are uniformly distributed in the tangential directions, so that the 
overall configuration is spherically—symmetric 10 . The corresponding shell action is [54] 

Sshell = - j rn eS dr , ml s = m 2 + L 2 /r 2 (r) , (3.2) 

where L is a parameter proportional to the angular momentum of the constituent particles. 
Its nonzero value provides a centrifugal barrier reflecting classical shells at low energies. 
Decreasing this parameter, we arrive to the regime of classical gravitational collapse. In 
what follows we switch between the scattering regimes by changing the parameter L = P” 1 . 
For completeness we derive the action (3.2) in Appendix A. 

Gravitational sector of the model is described by the Einstein-Hilbert action with the 
Gibbons-Hawking term, 

Seh = -T7T- [ d A Xy/^g R , (3.3) 

16vr J v 

S GH = ^~ [ k d 3 a y/\h\ (K - K 0 ) . (3.4) 

Here the metric g^ and curvature scalar R are defined inside the space-time volume V 
with the boundary 11 dV. The latter consists of a time-like surface at spatial infinity 
r = Too —> +oo and space-like surfaces at the initial and final times t = tij —> =poo. In 
Eq. (3.4) a are the coordinates on the boundary, h is the determinant of the induced metric, 
while K is the extrinsic curvature involving the outer normal. The parameter k equals +1 
(—1) at the time-like (space-like) portions of the boundary. To obtain zero gravitational 
action in flat space—time, we subtract the regulator Kq which is equal to the flat-space 
extrinsic curvature of the boundary [55]. For the sphere at infinity Kq = 2/roo, while the 
initial- and final-time hypersurfaces have Kq = 0. The Gibbons-Hawking term (3.4) will 
play an important role in our analysis. 

Let us first discuss the classical dynamics of the system. Equations of motion follow 
from variation of the total action 

S = S s heii + Seh + Sgh (3.5) 

with respect to the metric g^ u and the shell trajectory y a {r). In the regions inside and 
outside the shell the metric satisfies vacuum Einstein equations and therefore, due to 
Birkhoff theorem, is given by the flat and Schwarzschild solutions, respectively, see Fig. 3a. 
Introducing the spherical coordinates (t_, r ) inside the shell and Schwarzschild coordinates 
(t. |_, r) outside, one writes the inner and outer metrics in the universal form 

ds± = —f±(r)dt± + + r 2 dn 2 , (3.6) 

where 

/- = 1 , /+ = 1 - 2 M/r . (3.7) 

10 Similar models are used in astrophysics to describe structure formation [53]. 
lx We impose the Dirichlet boundary conditions on the variations of g^ at dV. 
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Figure 3. Gravitational field of the spherical dust shell: (a) in the original model, (b) in the model 
with modification at r « r w . 


The parameter M is the ADM mass which coincides with the total energy of the shell. In 
what follows we will also use the Schwarzschild radius I'h = 2 M. For the validity of the 
semiclassical approach we assume that the energy is higher than Planckian, M> 1. 

Equation for the shell trajectory is derived in Appendix B by matching the inner and 
outer metrics at the shell worldsheet with the Israel junction conditions [56, 57]. It can be 
cast into the form of an equation of motion for a particle with zero energy in an effective 
potential, 


r 2 + V eS (r) = 0 , (3.8) 

( 3 . 9 ) 

(L 2 + m 2 r 2 + 2 Mr 3 )" 

4r 4 (L 2 + m 2 r 2 ) 

This equation incorporates gravitational effects as well as the backreaction of the shell on 
the spacetime metric. The potential V e g(r) goes to —oo at r —> 0 and asymptotes to a 
negative value 12 1 — M 2 /m 2 at r = +oo, see Fig. 4. At large enough L the potential crosses 
zero at the points A and A' — the turning points of classical motion. A shell coming from 
infinity reflects from the point A back to r = +oo. When L decreases, the turning points 
approach each other and coalesce at a certain critical value 13 L = L*. At even smaller L 
the turning points migrate into the complex plane, see Fig. 5 (upper left panel), and the 
potential barrier disappears. Now a classical shell coming from infinity goes all the way to 
r = 0. This is the classical collapse. 

Now, we explicitly see an obstacle for finding the reflected semiclassical solutions at 
L < L* with the method of continuous deformations. Indeed, at large L the reflected 


12 Recall that the shell energy M is always larger than its rest mass m. 
12 For a massless shell (m = 0) the critical value is L * = 27A/ 2 /8. 
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A^eff r h 



Figure 4. Effective potential for the shell motion. 


solutions r = r(r) are implicitly defined as 


r(r) 


/ 


dr 


(3.10) 


V~V eS (r) 


where the square root is positive at r —> +00 + i 0. The indefinite integral is performed 
along the contour C running from r = +00 — iO to r = +00 + zO and encircling the turning 
point A — the branch point of the integrand (see the upper left panel of Fig. 5). As L is 
lowered, the branch point moves and the integration contour stays attached to it. However, 
at L = L* when the branch points A and A' coalesce, the contour C is undefined. It is 
therefore impossible to obtain reflected semiclassical solutions at L < L* from the classical 
solutions at L > L*. 

3.2 Modification 

To find physically relevant reflected trajectories at L < L*, we use the method of Sec. 2 and 
add an imaginary term itT %n i to the action. We consider Ti n t of the form (2.10), where the 
function w{r ) is concentrated in the vicinity of r = r w . The radius r w is chosen to be large 
enough, in particular, larger than the Schwarzschild radius r^ and the position r 2 4 of the 
right turning point A. Then the Einstein equations are modified only at r ~ r w , whereas 
the geometries inside and outside of this layer are given by the Schwarzschild solutions 
with masses M' and M, see Fig. 3b. To connect these masses, we solve the modified 
Einstein equations in the vicinity of r w . Inserting general spherically symmetric metric in 
the Schwarzschild frame, 


ds 2 = —f(t, r)dt 2 + -b r 2 dQ 2 , 


(3.11) 


into the ( tt ) component of Einstein equations, we obtain, 



(3.12) 
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e = 0: e > 0: 



Figure 5. Motion of the turning points in the complex r-plane. The points move along the arrows 
as L decreases. The cases of classical shell (upper left panel) and shell with the modified action 
(right panels) are shown. 

The solution reads 14 , 

2 M(r) ~ f +oa 

/ = 1-, where M(r) = M + ie / dr'w(r')F(f ) . (3.13) 

^ Jr 

This gives the relation 

M' = M + ie , e = e J dr' wF . (3-14) 

r^r w 

Here e > 0 is the new parameter of modification. As before, the ADM mass M of the 
system is conserved in the course of the evolution. It coincides with the initial and final 
energies of the shell which are, in turn, equal, as will be shown in Sec. 3.3, to the initial- and 
final-state energies in the quantum scattering problem. Thus, M is real, while the mass 
M' of the Schwarzschild space-time surrounding the shell acquires a positive imaginary 
part 1 ' 5 . The shell dynamics in this case is still described by Eq. (3.8), where M is replaced 
by M' in the potential (3.9). Below we find semiclassical solutions for small e > 0. In the 
end e will be sent to zero. 

Let us study the effect of the modification (3.14) on the semiclassical trajectories 
r = r(r) in Eq. (3.10). At L > L* the complex terms in V e s are negligible and the reflected 
trajectory is obtained with the same contour C as before, see the upper left panel of Fig. 5. 

14 The function / is time-independent due to the (tr) equation. 

15 In this setup the method of Sec. 2 is equivalent to analytic continuation of the scattering amplitude 
into the upper half-plane of complex ADM energy, cf. [28]. 
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The modification of V e g becomes important when L gets close to L* and the two turning 
points A and A' approach each other. Expanding the potential in the neighborhood of the 
maximum, we write, 

V eS (r) « V ma x - n 2 {r - r max ) 2 , (3.15) 

where V max , ^ and r max depend on L and M'. For real M' = M the extremal value 
Vmax is real and crosses zero when L crosses L*, whereas the parameters /t 2 > 0 and r max 
remain approximately constant. The shift of M' into the upper complex half-plane gives 
a negative imaginary part to V max , 

Irn Vmax = e + 0(e 2 ) < 0 , (3.16) 

where the last inequality follows from the explicit form (3.9). Now, it is straightforward to 
track the motion of the turning points using Eq. (3.15) as L decreases below L*. Namely, 
A and A' are shifted into the lower and upper half-planes as shown in Fig. 5 (upper right 
panel). Importantly, these points never coalesce. Physically relevant reflected solution at 
L < L* is obtained by continuously deforming the contour of integration in Eq. (3.10) while 
keeping it attached to the same turning point. As we anticipated in Sec. 2, a smooth branch 
of reflected semiclassical solutions parameterized by L exists in the modified system. 

If L is slightly smaller than L*, the relevant saddle-point trajectories reflect at 
Re rA > Th and hence never cross the horizon. A natural interpretation of the correspond¬ 
ing quantum transitions is over-barrier reflection from the centrifugal potential. However, 
as L decreases to L —> 0, the centrifugal potential vanishes. One expects that the semi¬ 
classical trajectories in this limit describe complete gravitational transitions proceeding via 
formation and decay of a black hole. 

After establishing the correspondence between the solutions at L > L* and L < L*, we 
take 16 e = 0. This leaves us with the original semiclassical equations which do not include 
the imaginary regularization terms. 

We numerically traced the motion of the turning point A as L decreases from large to 
small values, see Fig. 5 (lower panel). It approaches the singularity r = 0 at L — > 0. This 
behavior is confirmed analytically in Appendix C. Thus, at small L the contour C goes 
essentially along the real axis making only a tiny excursion into the complex plane near 
the singularity. It encircles the horizon r = from below. 

One remark is in order. For the validity of the low-energy gravity the trajectory should 
stay in the region of sub-Planckian curvature, R P i/\ p R p ‘ v> ' p ~ M 2 /r 6 <C 1. This translates 
into the requirement for the turning point 

M » M 1 / 3 . (3.17) 

On the other hand, we will see shortly that the dependence of the semiclassical amplitude 
on L drops off at ~ M 2 ; we always understand the limit L —> 0 in the sense of this 

inequality. Using the expressions for ?’n from Appendix C one verifies that the condition 
(3.17) can be satisfied simultaneously with in the semiclassical regime M> 1. 

lb This may be impossible in more complicated systems [27, 28, 30, 31] where the relevant saddle-point 
trajectories do not exist at 1 = 0. In that case one works at nonzero e till the end of the calculation. 
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Figure 6. The time contour corresponding to the semiclassical solution at small L. Solid and 
dashed lines correspond to interacting and free evolution respectively, cf. Fig. 2. 


3.3 5—matrix element 

The choice of the time contour. The action S reg entering the amplitude (1.3) is 
computed along the contour in complex plane of the asymptotic observer’s time t = t + . 
Since we have already found the physically relevant contour C for r(r), let us calculate the 
Schwarzschild time t+(r) along this contour. We write, 


f r dt + 
U(r ] = j dr — 


f r dr \/f+{r) - V e s(r) 

J V~ V es(r) f+(r ) 


(3.18) 


where the indefinite integral runs along C. In Eq. (3.18) we used the the definition of the 
proper time implying 


/+ 



1 1 
r 2 + /+ 


and expressed r 2 from Eq. (3.8). The integrand in Eq. (3.18) has a pole at the horizon 
r = r hi f+( r h ) = 0, which is encircled from below, see Fig. 5, lower panel. The half¬ 
residue at this pole contributes inr^ to t + each time the contour C passes close to it. The 
contributions have the same sign: although the contour C passes the horizon in the opposite 
directions, the square root in the integrand changes sign after encircling the turning point. 
Additional imaginary contribution comes from the integral between the real r-axis and the 
turning point A; this contribution vanishes at L —> 0. 

The image Ct of the contour C is shown in Fig. 6, solid line. Adding free evolution 
from t- 1_ = 0“ to t- 1_ = ti and from t+ = tf to t+ = 0 + (dashed lines), we obtain the contour 
analogous to the one in Fig. 2. One should not worry about the complex value of tf in 
Fig. 6: the limit tf +oo in the definition of S'—matrix implies that S reg does not depend 
on tf. Besides, the semiclassical solution r = r[t+) is an analytic function of t+ and the 
contour Ct can be deformed in complex plane as long as it does not cross the singularities 1 ' 
of r(t+). Below we calculate the action along Ct because the shell position and the metrics 
are real in the initial and final parts of this contour. This simplifies the calculation of the 
Gibbons-Hawking terms at t+ = ti and t+ = tf. 


17 In fact, Ct is separated from the real time axis by a singularity where r(t+) = 0. This is the usual 
situation for tunneling solutions in quantum mechanics and field theory [27, 28]. Thus, S reg cannot be 
computed along the contour in Fig. 2; rather, Ct or an equivalent contour should be used. 
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Interacting action. Now, we evaluate the action of the interacting system S(ti,tf ) 
entering S reg . We rewrite the shell action as 

Sshell — [ / T . n^eS > (3.19) 

JC V *eE 

where Eq. (3.10) was taken into account. The Einstein-Hilbert action (3.3) is simplified 
using the trace of the Einstein equations, R = — 87rT s/ieH ^, and the energy-momentum 
tensor of the shell computed in Appendix B, Eq. (B.5), 


Seh = 


rTf m 2 

dr - 

r 2rn e \\ 


f dr m 2 
lc V-Ves 2rn eS 


(3.20) 


An important contribution comes from the Gibbons-Hawking term at spatial infinity r = 
r oo —> +oo. The extrinsic curvature reads, 


K 


r oo 


rf+ + 4 /+ 

2 ry/U 


2 


Too 


# +0(r~ 3 ) 


(3.21) 


The first term here is canceled by the regulator Kq in Eq. (3.4). The remaining expression 
is finite at r^ -> +oo, 


Sgh 


m r m 

= 2 I dt+= 2 


dr V/+ " VeE 


V~Vefi f+ 


(3.22) 


where we transformed to integral running along the contour C using Eq. (3.18). Note that 
this contribution contains an imaginary part 


Im S G h 


M t . 

= -y Im(*/ - U 


(3.23) 


Finally, in Appendix D we evaluate the Gibbons-Hawking terms at the initial- and final¬ 
time hypersurfaces. The result is 


Sgh 




VM 2 - m 2 ^ ( M (2M 2 - m 2 ) 

2 + 4VM 2 - m 2 " 


(3.24) 


where r,j are the radii of the shell at the endpoints of the contour C. The latter radii are 
real, and so are the terms (3.24). 

Summing up the above contributions, one obtains, 


S(ti,t f ) 


f dr 

m 2 — 2 m 2 s 

M a//+ — V e ff 

1 V-VeS 

2m e ff 

2 U _ 


c 


+ 


V M 2 — m 2 , 

- 5 - (n + r f ) + 


M (2M 2 —m 2 ) 
2 yfM 2 -m 2 


(3.25) 


This expression contains linear and logarithmic divergences when r, j are sent to infinity. 
Note that the divergences appear only in the real part of the action and thus affect only 
the phase of the reflection amplitude but not its absolute value. 
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Initial and final—state contributions. The linear divergence in Eq. (3.25) is related 
to free motion of the shell in the asymptotic region r —> +oo, whereas the logarithmic one 
is due to the 1/r tails of the gravitational interaction in this region. Though the 1/r terms 
in the Lagrangian represent vanishingly small gravitational forces in the initial and final 
states, they produce logarithmic divergences in S(ti,tj ) when integrated over the shell 
trajectory. To obtain a finite matrix element, we include 18 these terms in the definition of 
the free action Sq. In Appendix E the latter action is computed for the shell with energy 
M, 

i-n rri 

So(0~,ti)= / pi(r, M)dr - Mti , S 0 (t fl O + )= / p f (r,M)dr + Mtf , (3.26) 

J n Jrf 


where n .2 are the positions of the shell at t+ = O 18 and 


Pi,f(r,M) = =f 


\J M 2 — m 2 + 


M(2M 2 - m 2 ) 
2 r\J M 2 — m 2 


(3.27) 


are the initial and final shell momenta with 1/r corrections. 

The path integral (2.2) for the amplitude involves free wavefunctions Tj(ri) and Ty/^) 
of the initial and final states. We consider the semiclassical wavefunctions of the shell with 
fixed energy E, 




'lq(ri) ~ exp < i / pi(r',E)dr' 


' r o 


rr 2 


— exp i I pf(r',E)dr' 
I r 0 


(3.28) 


where p-j j are the same as in Eq. (3.27). In fact, the energy E is equal to the energy of the 
semiclassical solution, E = M. Indeed, the path integral (2.2) includes integration over 
the initial and final configurations of the system, i.e. over r\ and r 2 in the shell model. The 
condition for the stationary value of r\ reads, 

d d 

— logTj + i—S 0 (0~,ti) = 0 =» pi(ri,E) = pi(n,M) , (3.29) 

or i ari 

and similarly for r 2 - This implies equality of E and M. Note that the parameter ro 2 M 
in Eq. (3.28) fixes the phases of the initial and final wavefunctions. Namely, the phases 
vanish at r = tq. 


The result. Combining the contributions (3.26), (3.28) with Eq. (3.25), we obtain the 
exponent of the 5-matrix element (2.2), 


Sreg — 


dr 


m 2 - 2?B 2 ff M y/f+ - Vetf 


J V~VeS [ 2 r n eS 2 
c 

- y/M 2 - m 2 r r i± r l _ 2r 0 ^ + 


U 

M(2M 2 - m 2 ) 


1 - log(rjrj/rg)^ 


(3.30) 


2 y/M 2 - m 2 

Notice that this result does not depend on any regularization parameters. It is straightfor¬ 
ward to check that Eq. (3.30) is finite in the limit r,j —> +oo. In Fig. 7 we plot its real and 


s This relies on the freedom in splitting the total Lagrangian into “free” and “interacting” parts. 
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Figure 7. Real (a) and imaginary (b) parts of S reg at m = 0 as functions of L. For the real part 
we take ?'o = 10 M. 


imaginary parts as functions of L for the case of massless shell (m = 0). The imaginary 
part vanishes for the values L > L* corresponding to classical reflection. At smaller L the 
imaginary part is positive implying that the reflection probability 

'Pfi — g — 2I m >S reg 


is exponentially suppressed. Importantly, S reg does not receive large contributions from 
the small-r region near the spacetime singularity. It is therefore not sensitive to the effects 
of trans-Planckian physics. 

In the most interesting case of vanishing centrifugal barrier L —> 0 the only imaginary 
contribution to S reg comes from the residue at the horizon = 2 M in Eq. (3.30), recall 
the contour C in Fig. 5. The respective value of the suppression exponent is 


21m Sj'gg 


= 27rMRes r=r . 

L =0 


VU - VeS 

V=v*f+ 


nrl • 


(3.31) 


This result has important physical implications. First, Eq. (3.31) depends only on the total 
energy M of the shell but not on its rest mass m. Second, the suppression coincides with 
the Bekenstein-Hawking entropy of a black hole with mass M. The same suppression was 
obtained in [39, 40] for the probability of emitting the total black hole mass in the form of a 
single shell. We conclude that Eq. (3.31) admits physical interpretation as the probability 
of the two-stage reflection process where the black hole is formed in classical collapse 
with probability of order 1, and decays afterwards into a single shell with exponentially 
suppressed probability. 

One may be puzzled by the fact that, according to Eq. (3.30), the suppression receives 
equal contributions from the two parts of the shell trajectory crossing the horizon in the 
inward and outward directions. Note, however, that the respective parts of the integral 
(3.30) do not have individual physical meaning. Indeed, we reduced the original two- 
dimensional integral for the action to the form (3.30) by integrating over sections of constant 
Schwarzschild time. Another choice of the sections would lead to an expression with a 
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different integrand. In particular, using constant-time slices in Painleve or Finkelstein 
coordinates one obtains no imaginary contribution to S reg from the inward motion of the 
shell, whereas the contribution from the outward motion is doubled. The net result for the 
probability is, of course, the same. 19 

The above result unambiguously shows that the shell model, if taken seriously as a 
full quantum theory, suffers from the information paradox. Indeed, transition between 
the only two asymptotic states in this theory — contracting and expanding shell — is 
exponentially suppressed. Either the theory is intrinsically non-unitary or one has to take 
into consideration an additional asymptotic state of non-evaporating eternal black hole 
formed in the scattering process with probability 1 — Vf t . 

On the other hand, the origin of the exponential suppression is clear if one adopts 
a modest interpretation of the shell model as describing scattering between the narrow 
wavepackets in field theory. Hawking effect implies that the black hole decays predomi¬ 
nantly into configurations with high multiplicity of soft quanta. Its decay into a single hard 
wavepacket is entropically suppressed. One can therefore argue [42] that the suppression 
(3.31) is compatible with unitarity of field theory. However, the analysis of this section is 
clearly insufficient to make any conclusive statements in the field theoretic context. 

As a final remark, let us emphasize that besides the reflection probability our method 
allows one to calculate the phase of the scattering amplitude Re S reg . At L = m = 0 it can 
be found analytically, 

Re S re g = 2Mr 0 + 2M 2 log(r 0 /2 M) + M 2 . (3.32) 

Notice that apart from the trivial first term, the phase shift is proportional to the entropy 
B oc Af 2 ; this is compatible with the dependence conjectured in [20]. The phase (3.32) 
explicitly depends on the parameter vq of the initial- and final-state wavefunctions. 

3.4 Relation to the Hawking radiation 

In this section we deviate from the main line of the paper which studies transitions between 
free-particle initial and final states, and consider scattering of a shell off an eternal pre¬ 
existing black hole. This will allow us to establish a closer relation of our approach to the 
results of [39, 40] and the Hawking radiation. We will focus on the scattering probability 
and thus consider only the imaginary part of the action. 

The analysis essentially repeats that of the previous sections with several differences. 
First of all, the inner and outer space-times of the shell are now Schwarzschild with the 
metric functions 


/_ = 1 - 2 M BH /r , /+ = 1 - 2 (M bh + M)/r , (3.33) 

where Mbh is the eternal black hole mass and M denotes, as before, the energy of the 
shell. The inner and outer metrics possess horizons at r7 = 2Mbh and = 2(A Ibh + M), 
respectively. The shell motion is still described by Eq. (3.8), where the effective potential is 
obtained by substituting expressions (3.33) into the first line of Eq. (3.9). Next, the global 

19 Note that our semiclassical method is free of uncertainties [58-60] appearing in the approach of [39]. 
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Figure 8. Penrose diagram for scattering of a classical shell from an eternal black hole. 


space-time has an additional boundary r = r'^ -» +oo at the second spatial infinity of 
the eternal black hole, see Fig. 8. We have to include the corresponding Gibbons-Hawking 
term, cf. Eq. (3.22), 


Sgh 



(3.34) 


Finally, the eternal black hole in the initial and final states contributes into the free action 
So- We use the Hamiltonian action of an isolated featureless black hole in empty space- 
time [61], 

Sq,bh = -M B h J dt+ , (3.35) 


where, as usual, the time variable coincides with the asymptotic time 20 . Since we do not 
equip the black hole with any degrees of freedom 21 , its initial- and final-state wavefunctions 
are \k bh = 1- 

Adding new terms (3.34), (3.35) to the action (3.30), one obtains 22 , 


Im S re g 



dr 


m 2 - 2m 2 ff M B h + M a//+ - V e s 
2 m eS 2 /+ 


Mbh a //- ~ VsS 
2 /- 


(3.36) 

where the integration contour C is similar to that in Fig. 5 (lower panel), it bypasses the 
two horizons r7 and in the lower half of complex r-plane. In the interesting limit of 
vanishing centrifugal barrier L —> 0 the imaginary part of the action is again given by the 
residues at the horizons, 


21m S, 


reg 


L^2n(M BH + M) 


= [( r t ) 2 - (r h ) 2 } 


27 tM B h Res - 

h 


Vf- - VeS 


(3.37) 


20 To compute the phase of the scattering amplitude, one should also take into account long-range inter¬ 
action of the black hole with the shell in the initial and final states. 

' 21 The choice of a model for an isolated black hole is the main source of uncertainties in scattering problems 
with eternal black holes. 

22 Recall that the free action enters this formula with the negative sign, So,bh = — ti). Note also 

that the original Gibbons-Hawking term (3.22) is proportional to the total ADM mass Mbh + M. 
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Interpretation of this result is similar to that in the previous section. At the first stage of 
transition the black hole swallows the shell with probability of order one and grows to the 
mass Mbh + M. Subsequent emission of the shell with mass M involves suppression 

V fi e~ B + +B ~ , (3.38) 

where B± = are the entropies of the intermediate and final black holes. This 

suppression coincides with the results of [39, 40]. 

At Mbh = 0 the process of this section reduces to reflection of a single self-gravitating 
shell and expression (3.37) coincides with Eq. (3.31). In the other limiting case M <C Mbh 
the shell moves in external black hole metric without back reaction. Reflection probability 
in this case reduces to the Boltzmann exponent 

V fi - e~ M ' TH , 

where we introduced the Hawking temperature Th = 1/(87 tMbh)- One concludes that 
reflection of low-energy shells proceeds via infall into the black hole and Hawking evapo¬ 
ration, whereas at larger M the probability (3.38) includes back-reaction effects. 

3.5 Space—time picture 

Let us return to the model with a single shell considered in Secs. 3.1-3.3. In the previous 
analysis we integrated out the non-dynamical metric degrees of freedom and worked with 
the semiclassical shell trajectory (t + (r), r(r)). It is instructive to visualize this trajectory in 
regular coordinates of the outer space-time. Below we consider the case of ultrarelativistic 
shell with small angular momentum: L —» 0 and M m. One introduces Kruskal 
coordinates for the outer metric, 

U = —(r/2M' - 1) 1/2 e ( r - f +)/ 4M ' s v = (r/2M / - 1) 1/2 e (r+t +)/ 4M ' . (3.39) 

We choose the branch of the square root in these expressions by recalling that M' differs 
from the physical energy M by an infinitesimal imaginary shift, see Eq. (3.14). The initial 
part of the shell trajectory from t + = ti to the turning point A (Figs. 5, 6) is approximately 
mapped to a light ray V = Vo > 0 as shown in Fig. 9. Note that in the limit L —> 0 the 
turning point A is close to the singularity r = 0, but does not coincide with it. At the 
turning point the shell reflects and its radius r(r) starts increasing with the proper time r. 
This means that the shell now moves along the light ray U = Uq > 0, and the direction of r 
is opposite to that of the Kruskal time U+V. The corresponding evolution is represented by 
the interval ( A,tf) in Fig. 9. We conclude that at t+ = tf the shell emerges in the opposite 
asymptotic region in the Kruskal extension of the black hole geometry. This conclusion may 
seem puzzling. However, the puzzle is resolved by the observation that the two asymptotic 
regions are related by analytic continuation in time. Indeed it is clear from Eqs. (3.39) that 
the shift £_|_ i —> t+ — 4irMi corresponds to total reflection of Kruskal coordinates U —> —U, 
V — > — V. Precisely this time—shift appears if we extend the evolution of the shell to the 
real time axis (point t'j in Fig. 6). At t+ = t/j the shell emerges in the right asymptotic 
region 23 with future-directed proper time r. The process in Fig. 9 can be viewed as a 

Although the shell coordinate r(t/) is complex, we identify the asymptotic region using Rer(t/) which 
is much larger than the imaginary part. 
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Figure 9. Trajectory of the shell in Kruskal coordinates of the outer metric. Black dashed lines 
show the singularities UV = 1. 


shell-antishell annihilation which is turned by the analytic continuation into the transition 
of a single shell from fj to t'^. 

Now, we write down the space-time metric for the saddle-point solution at m = 0 
and L —> 0. Recall that in this case the shell moves along the real r-axis. We therefore 
introduce global complex coordinates (r, t + ), where t + belongs to Ct and r is real positive. 
The metric is given by analytic continuation of Eqs. (3.6), (3.7), 


ds 2 


~ f 1 “ idSfcl) " dt l + dr 2 + rW ’ r < r shell(t+) 

~ (1 ~ dt\ + + r 2 dQ 2 , r > r she u(t+) , 

\ / r 


(3.40) 


where we changed the inner time f_ to by matching them at the shell worldsheet 
r = r she ii(t + ). Importantly, the metric (3.40) is regular at the origin r = 0 which is never 
reached by the shell. It is also well defined at = 2 M due to the imaginary part of 
M'\ in the vicinity of the Schwarzschild horizon r/j the metric components are essentially 
complex. Discontinuity of Eq. (3.40) at r = r s h e ii(t +) is a consequence of the d-function 
singularity in the shell energy-momentum tensor. This makes the analytic continuation of 
the metric ill-defined in the vicinity of the shell trajectory. We expect that this drawback 
disappears in the realistic field-theory setup where the saddle-point metric will be smooth 
(and complex-valued) in Schwarzschild coordinates. 


4 Massless shell in AdS 
4.1 Reflection probability 

In this and subsequent sections we subject our method to further tests in more complicated 
shell models. Here we consider a massless shell in 4-dimensional AdS space-time. The 
analysis is similar to that of Sec. 3, so we will go fast over details. 

The shell action is still given by Eq. (3.2) with m e ff = T/r, while the Einstein-Hilbert 
action is supplemented by the cosmological constant term, 

Seh = / d 4 xy/^g (R ~ 2A) . (4.1) 

167T J v 
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Here A = —3 /l 2 , l is the AdS radius. The Gibbons-Hawking term has the form (3.4), 
where now the regulator at the distant sphere 


K n 


— t H— 


(4.2) 


is chosen to cancel the gravitational action of an empty AdS,*. 
outside the shell is AdS and AdS-Schwarzschild, respectively, 


/- - 1 + ^ , 


2 M r 2 

f+ ~ 1-b V7 

r R 


The metric inside and 


(4.3) 


where M is the shell energy. The trajectory of the shell obeys Eq. (3.8) with the effective 
potential given by the first line of Eq. (3.9), 


= 1 T ^ 


(L 2 + 2 Mr 3 ) 2 
4r 4 L 2 


(4.4) 


The e-modification again promotes M in this expression to M' = M + ie. Repeating the 
procedure of Sec. 3.2, we start from the reflected trajectory at large L. Keeping e > 0, we 
trace the motion of the turning point as L decreases 21 . The result is a family of contours 
C spanned by the trajectory in the complex r-plane. These are similar to the contours in 
Fig. 5. In particular, at L —> 0 the contour C mostly runs along the real axis encircling the 
AdS-Schwarzschild horizon from below, as in the lower panel of Fig. 5. 

Calculation of the action is somewhat different from that in flat space. First, the 
space-time curvature is now non-zero everywhere. Trace of the Einstein’s equations gives 25 
R = 4A. The Einstein-Hilbert action takes the form, 


Seh 


A 

2 


r shell 



0 


T 00 



T shell 


= J Jl 2 ( dt + ~ dt ~) 

shell 



dt+ . (4.5) 


The last term diverging at —> 00 is canceled by the similar contribution in the Gibbons- 

Hawking term at spatial infinity, 

Sgh I- = (w~ y ) i dt+ ■ (46) 

Second, unlike the case of asymptotically flat space-time, Gibbons-Hawking terms at the 
initial- and final-time hypersurfaces t+ = tij vanish, see Appendix D. Finally, the canon¬ 
ical momenta 26 of the free shell in AdS, 

M 2 L 2 

Plf(r, M) = (1 + r 2 /p) 2 " r 2(l +r 2/ Z 2) ’ ( 4J ) 

are negligible in the asymptotic region r —> + 00 . Thus, the terms involving pij in the free 
action (3.26) and in the initial and final wavefunctions (3.28) are vanishingly small if the 

24 Alternatively, one can start from the flat-space trajectory and continuously deform it by introducing 
the AdS radius l. 

25 In the massless case the trace of the shell energy-momentum tensor vanishes, T shM = 0, see Eqs. (B.3). 
26 They follow from the dispersion relation m 2 ff = — g 00 M 2 — g rr p 2 . 
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normalization point ?’o is large enough. This leaves only the temporal contributions in the 
free actions, 


So(0 ,U) + S 0 (tf,0 + ) = M(tf - U) . 


(4.8) 


Summing up Eqs. (4.5), (4.6), (4.8) and the shell action (3.2), we obtain, 


Sreg — 


dr 


IcV^v^lu 2 


Vf+ - Kff Vf-~ V eA + M V / / + -Eeff 


/+ 


f- 


f+ 


L 

r 


^ 2 L ‘ t) + 5/d at i ^°’ (4 ' 9) 

where the integration contour in the last expression goes below the pole at r = . The 

integral (4.9) converges at infinity due to fast growth of functions /+ and /_. In particular, 
this convergence implies that there are no gravitational self-interactions of the shell in the 
initial and final states due to screening of infrared effects in AdS. 

The imaginary part of Eq. (4.9) gives the exponent of the reflection probability. It is 
related to the residue of the integrand at r^, 


2Inr S reg = 27r ( =) + M 


l 2 


1 


BXSr - rk Mr) 


= 


(4.10) 


We again find that the probability is exponentially suppressed by the black hole entropy. 
Remarkably, the dependence of the reflection probability on the model parameters has 
combined into which is a complicated function of the AdS-Schwarzschild parameters M 
and l. 


4.2 AdS/CFT interpretation 

Exponential suppression of the shell reflection has a natural interpretation within the 
AdS/CFT correspondence [9, 62, 63]. The latter establishes relationship between gravity 
in AdS and strongly interacting conformal held theory (CFT). Consider three-dimensional 
CFT on a manifold with topology Mx S 2 parameterized by time t and spherical angles 6. 
This is the topology of the AdS 4 boundary, so one can think of the CFT 3 as living on this 
boundary. Let us build the CFT dual for transitions of a gravitating shell in AdS 4 . Assume 
the CFT 3 has a marginal scalar operator 0(t, 0); its conformal dimension is A = 3. This 
operator is dual to a massless scalar held cj) in AdS 4 . 

Consider now the composite operator 

Om(^o) = exp | J dtd 2 6G M {t-to)d(t,0)\ , (4.11) 

where Gm{P) is a top-hat function of width At 3> 1/M. This operator is dual to a spherical 
wavepacket (coherent state) of the <)>-held emitted at time to from the boundary towards 
the center of AdS [64, 65]. With an appropriate normalization of Gm {£)> the energy of the 
wavepacket is M. Similarly, the operator 0/^(to) is dual to the wavepacket absorbed on 
the boundary at time to- Then the correlator 2, 

GM = (d+(ri)d M m (4.12) 

2 'The time nl is needed for the wavepacket to reach the center of AdS and come back if it moves 
classically [65]. 
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Figure 10. Conformal diagram for scattering of a massless shell in AdS 4 . Creation and annihilation 
of the shell at the AdS boundary correspond to insertions of composite operators in the CFT dual. 


is proportional to the amplitude for reflection of the contracting wavepacket back to the 
boundary. If the width of the wavepacket is small enough, At <C l, the 0—field can be treated 
in the eikonal approximation and the wavepacket follows a sharply defined trajectory. In 
this way we arrive to the transition of a massless spherical shell in AdS 4 , see Fig. 10. 

Exponential suppression of the transition probability implies respective suppression of 
the correlator (4.12). However, the latter suppression is natural in CFT 3 because the state 
created by the composite operator Om{ 0) is very special. Submitted to time evolution, 
it evolves into a thermal equilibrium which poorly correlates with the state destroyed by 
0 +(ttZ). Restriction of the full quantum theory in AdS .4 to a single shell is equivalent 
to a brute-force amputation of states with many soft quanta in unitary CFT 3 . Since the 
latter are mainly produced during thermalization, the amputation procedure leaves us with 
exponentially suppressed 5-matrix elements. 

5 Charged shells 
5.1 Elementary shell 

Another interesting extension of the shell model is obtained by endowing the shell with 
electric charge. The corresponding action is the sum of Eq. (3.5) and the electromagnetic 
contribution 

Sem = -T 7 T- / d 4 x^g Ff w -Q A a dy a , (5.1) 

!6vr J ^ J she u 

where A M is the electromagnetic field with stress tensor F /1V = d^Ay — dyA^ and Q is the 
shell charge. This leads to Reissner-Nordstrom (RN) metric outside the shell and empty 
flat space-time inside, 

/+ = — + A 0+ = — ; /_ = 1, A 0 _=0. (5.2) 

rjr* iyt Zj rj™ 

Other components of A^ are zero everywhere. Importantly, the outside metric has two 
horizons 

rj l ±) = M ± \JM' 2 - Q 2 (5.3) 
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Figure 11. Motion of the turning points and the contour C defining the trajectory for (a) the 
model with elementary charged shell and (b) the model with discharge. 


at Q < M. At Q > M the horizons lie in the complex plane, and the shell reflects 
classically. Since the latter classical reflections proceed without any centrifugal barrier, we 
set L = 0 henceforth. The semiclassical trajectories will be obtained by continuous change 
of the shell charge Q. 

The evolution of the shell is still described by Eq. (3.8) with the effective potential 
constructed from the metric functions (5.2), 


Kff = 1 - 


(m 2 — Q 2 + 2 Mr) 2 
4 m, 2 r 2 


This potential always has two turning points on the real axis, 


E4,A' 


Q 2 — m 2 
2(M =F m) 


(5.4) 


(5.5) 


The shell reflects classically from the rightmost turning point at Q > M. In the opposite 
case Q < M the turning points are covered by the horizons, and the real classical solutions 
describe black hole formation. 

We find the relevant semiclassical solutions at Q < M using e-modification. Since 
the modification term (2.10) does not involve the electromagnetic field, it does not affect 
the charge Q giving, as before, an imaginary shift to the mass, + ie. A notable 

difference from the case of Sec. 3 is that the turning points (5.5) are almost real at Q < M. 
The semiclassical trajectories therefore run close to the real r-axis 28 for any Q. On the 
other hand, the horizons (5.3) approach the real axis with Imr[ + ^ > 0 and Imr[ ^ < 0 as Q 
decreases. Thus, the saddle-point trajectories are defined along the contour C in Fig. 11a 
bypassing r h + ' and r j, ^ from below and from above, respectively. 

Since the semiclassical motion of the shell at Q < M proceeds with almost real r(r), 
we can visualize its trajectory in the extended RN geometry, see Fig. 12. The shell starts 
in the asymptotic region X, crosses the outer and inner horizons and \ repels from 
the time-like singularity due to electromagnetic interaction, and finally re-emerges in the 
asymptotic region X'. At first glance, this trajectory has different topology as compared 
to the classical reflected solutions at Q > M: the latter stay in the region X at the final 

28 The overall trajectory is nevertheless complex because t+ £ C, see below. 
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Figure 12. Conformal diagram for the extended Reissner-Nordstrom space-time. Semiclassical 
trajectory running along the contour C in Fig. 11a is shown by the red line. The grey region does 
not exist in theories with dynamical charged fields. 


time t+ = tf. However, following Sec. 3.5 we recall that the Schwarzschild time t+ of the 
semiclassical trajectory is complex in the region T ', 


Im (tf - ti) 


2 vr 

/;<>i +) > 


2 vr 


/+( r 



(5.6) 


where we used Eq. (3.18) and denoted by ti and tf the values of t + at the initial and 
final endpoints of the contour C in Fig. 11a. Continuing tf to real values, we obtain the 
semiclassical trajectory arriving to the region Z in the infinite future 29 , cf. Sec. 3.5. This 
is what one expects since the asymptotic behavior of the semiclassical trajectories is not 
changed in the course of continuous deformations. 

Let us now evaluate the reflection probability. Although the contour C is real, it 
receives imaginary contributions from the residues at the horizons. Imaginary part of the 
total action comes 30 from Eq. (3.30) and the electromagnetic term (5.1). The latter takes 
the form, 

Sem = ~ J d A x^g(- ^ J d A Xy/^gF 2 u , (5.7) 

where we introduced the shell current j 11 . used Maxwell equations = 47 rj^ and 

integrated by parts. From Eq. (5.2) we find, 


Sem = ^ 


oo 

I dt+ / 


r 2 dr( - 


Tshell 


Q 2 f dt + 

~2 J — 

shell 


(5.8) 


29 Indeed, the coordinate systems that are regular at the horizons and are periodic in the 

imaginary part of t+ with periods 47ri//j_(r( + ') and 47ri//j_(r^ _> ), respectively. Analytic continuation to 
real tf implies shifts by half-periods in both systems, see Eq. (5.6). This corresponds to reflections of the 
trajectory final point in Fig. 12 with respect to points O and 0^~\ 

30 Gibbons-Hawking terms at t = Uj are different in the case of charged shell from those in Eq. (3.30). 
However, they are real and do not contribute into ImStot- 
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Combining this with Eq. (3.30), we obtain, 


2ImS' re5 =Im 


dr 


'C V~VeK . 


— m + 


M _OlWU^VeZ 


=2tt [Res (+) - Res r=r (-) 


/+ 

r ) /+(r) 


= 7r ( r i +) ) 2 - 7r ( r ) l ’) 


H\2 


(5.9) 


After non-trivial cancellation we again arrive to a rather simple expression. However, this 
time 21m Stot is not equal to the entropy of the RN black hole, -Brat = 7 r ( r i + J- 

The physical interpretation of this result is unclear. We believe that it is an artifact 
of viewing the charged shell as an elementary object. Indeed, in quantum mechanics of an 
elementary shell the reflection probability should vanish at the brink Q = M of classically 
allowed transitions. It cannot be equal to Brat which does not have this property unlike 
the expression (5.9). We now explain how the result is altered in a more realistic setup. 


5.2 Model with discharge 

Recall that the inner structure of charged black holes in theories with dynamical fields is 
different from the maximal extension of the RN metric. Namely, the RN Cauchy horizon 
^ suffers from instability due to mass inflation and turns into a singularity [44-47]. 
Besides, pair creation of charged particles forces the singularity to discharge [43, 48, 49]. 
As a result, the geometry near the singularity resembles that of a Schwarzschild black hole, 
and the singularity itself is space-like. The part of the maximally extended RN space-time 
including the Cauchy horizon and beyond (the grey region in Fig. 12) is never formed in 
classical collapse. 

Let us mimic the above discharge phenomenon in the model of a single shell. Although 
gauge invariance forbids non-conservation of the shell charge Q, we can achieve essentially 
the same effect on the space-time geometry by switching off the electromagnetic interaction 
at r —> 0. To this end we assume spherical symmetry and introduce a dependence of the 
electromagnetic coupling on the radius 31 . This leads to the action 

= hS iixV ^~ Q S A ‘ iy ° ’ <5 - 10) 

shell 

where e(x) is a positive form-factor starting from e = 0 at x = 0 and approaching e —> 1 
at x —> Too. We further assume 

e(x) < x , (5-11) 

the meaning of this assumption will become clear shortly. Note that the action (5.10) is 
invariant under gauge transformations, as well as diffeomorphisms preserving the spherical 
symmetry. The width of the form-factor e(r/Q) in Eq. (5.10) scales linearly with Q to 
mimic larger discharge regions at larger Q. 

The new action (5.10) leads to the following solution outside the shell, 

OO 

2 MO /' p 2 (r') 

/+ = 1- 1 - air IQ) , Ao+ = a{r/Q) , where a(x) = / -T- dx f . (5.12) 

r r J xJ 

X 

31 Alternatively, the discharge can be modeled by introducing nonlinear dielectric permittivity [66]. 
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The space-time inside the shell is still empty and flat. As expected, the function /+ 
corresponds to the RN metric at large r and the Schwarzschild one at r —> 0. Moreover, 
the horizon r^ satisfying /+(?’/ r ) = 0 is unique due to the condition (5.11). It starts from 
= 2M at Q = 0, monotonically decreases with Q and reaches zero at Q* = 2M/o(0). 
At Q > Q* the horizon is absent and the shell reflects classically. 

The subsequent analysis proceeds along the lines of Secs. 3, 4. One introduces effective 
potential for the shell motion, cf. Eq. (5.4), 


= 1 - 


(2 Mr — Qra(r/Q) + m 2 ) 2 
4 m 2 r 2 


(5.13) 


introduces e-regularization, M h > M' = M + ie, and studies motion of the turning points 
of the shell trajectory as Q decreases. If M, Q m, this analysis can be performed for 
general e(x). In this case the relevant turning point va is real and positive for Q > Q*. 
At Q ~ Q* it comes to the vicinity of the origin r = 0 where the function a{r/Q ) can be 
expanded up to the linear term. The position of the turning point is 



M + m + 


a(0)Q 

2 


+ 



(5.14) 


where b 2 = 
interval 


—da/dx | t=q is positive according to Eq. (5.12). As Q decreases within the 


Q* 


2m(l — b) 

«( 0 ) 


> Q > Q* 


2m(l + b) 

«( 0 ) 


(5.15) 


the turning point makes an excursion into the lower half of the r-plane, goes below the 
origin and returns to the real axis on the negative side, see Fig lib. For smaller charges 
ta is small and stays on the negative real axis. The contour C defining the trajectory is 
shown in Fig. lib. It bypasses the horizon r from below, goes close to the singularity, 
encircles the turning point and returns back to infinity. This behavior is analogous to that 
in the case of neutral shell. 

Finally, we evaluate the imaginary part of the action. The electromagnetic contribution 
is similar to Eq. (5.8), 



(5.16) 


shell 


However, in contrast to Sec. 5.1, the trace of the gauge field energy-momentum tensor does 
not vanish due to explicit dependence of the gauge coupling on r (cf. Eq. (B.3b)), 


rpt n _ de 
EM n 87re 3 dr 


(5.17) 


This produces non-zero scalar curvature R = —8ttT em ^ in the outer region of the shell, 
and the Einstein-Hilbert action receives an additional contribution, 


ASeh 



r de 
e 3 dr 



shell 


dt+ 


Q 2 e 2 

4 r 



(5.18) 
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where in the second equality we integrated by parts. Combining everything together, we 
obtain (cf. Eq. (5.9)), 


21m S reg = Im 


dr 

— m + 

Jm- 

Q 2 e 

V-v eS 

2 r 

r=r h 

Q 2 e 2 

Qa 

\ 1 

2 r 

2~ y 

'/+ 


Qa \ VJ- 


-v eS 


U 


= nr 


h » 


(5.19) 


where non-trivial cancellation happens in the last equality for any e(x). To sum up, we 
accounted for the discharge of the black hole singularity and recovered the intuitive result: 
the reflection probability is suppressed by the entropy of the intermediate black hole 32 . 


6 Conclusions and outlook 

In this paper we proposed a semiclassical method to calculate the 5-matrix elements for the 
two-stage transitions involving collapse of multiparticle states into a black hole and decay of 
the latter into free particles. Our semiclassical approach does not require full quantization 
of gravity. Nevertheless, it consistently incorporates backreaction of the collapsing and 
emitted quanta on the geometry. It reduces evaluation of the 5-matrix elements to finding 
complex-valued solutions of the coupled classical Einstein and matter field equations with 
certain boundary conditions. 

An important technical ingredient of the method is the regularization enforcing the 
semiclassical solutions to interpolate between the in- and out- asymptotic states consist¬ 
ing of free particles in flat space-time. As a consequence, one works with the complete 
semiclassical solutions describing formation and decay of the intermediate black hole. This 
distinguishes our approach from the standard semiclassical expansion in the black hole 
background. In addition, the same regularization allows us to select the relevant semi¬ 
classical configurations by continuous deformation of the real solutions describing classical 
scattering at lower energies. The final result for the 5-matrix elements does not dependent 
on the details of the regularization. 

We illustrated the method in a number of toy models with matter in the form of thin 
shells. We have found that the relevant semiclassical solutions are complex-valued and 
defined in the complexified space-time in the case of black hole-mediated processes. They 
avoid the high-curvature region near the black hole singularity thus justifying our use of 
the semiclassical low-energy gravity. In particular, the Planck-scale physics near the black 
hole singularity is irrelevant for the processes considered in this paper. 

The method has yielded sensible results for transition amplitudes in the shell mod¬ 
els. Namely, we have found that the probabilities of the two-stage shell transitions are 
exponentially suppressed by the Bekenstein-Hawking entropies of the intermediate black 
holes. If the shell model is taken seriously as a full quantum theory, this result implies that 
its 5-matrix is non-unitary. However, the same result is natural and consistent with uni- 
tarity if the shells are interpreted as describing scatterings of narrow wavepackets in field 

i2 We do not discuss the phase of the scattering amplitude as it essentially depends on our choice of the 
discharge model. 
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theory. The exponential suppression appears because we consider a very special exclusive 
process: formation of a black hole by a sharp wavepacket followed by its decay into the 
same packet. Our result coincides with the probability of black hole decay into a single shell 
found within the tunneling approach to Hawking radiation [39, 40] and is consistent with 
interpretation of the Bekenstein-Hawking entropy as the number of black hole microstates 
[42] . Considering the shell in AdS 4 space-time we discussed the result from the AdS/CFT 
viewpoint. 

In the case of charged shells our method reproduces the entropy suppression only if 
instability of the Reissner-Nordstrom Cauchy horizon with respect to pair-production of 
charged particles is taken into account. This suggests that the latter process is crucial for 
unitarity of transitions with charged black holes at the intermediate stages. 

Besides the overall probability, our method yields the phase of the transition amplitude. 
The latter carries important information about the scattering process, in particular, about 
its initial and final states. In the case of a neutral shell in asymptotically flat space-time the 
phase contains a logarithmically divergent term due to long-range Newtonian interactions 
and terms proportional to the black hole entropy. This is consistent with the behavior 
conjectured in [20]. 

We consider the above successes as an encouraging confirmation of the viability of our 
approach. 

The shell models are too simple to address many interesting questions about the black 
hole 5-matrix. These include the expected growth of the transition probability when the 
final state approaches the Hawking-like state with many quanta, and the recent conjecture 
about sensitivity of the amplitudes to small changes in the initial and final states [67]. A 
study of these issues will require application of our method to a genuinely field-theoretic 
setup. Let us anticipate the scheme of such analysis. As an example, consider a scalar field 
4> minimally coupled to gravity. For simplicity, one can restrict to transitions between initial 
and final states with particles in the s-wave. These states are invariant under rotations. 
Then the respective semiclassical solutions are also expected to possess this symmetry. 
They satisfy the complexified wave- and Einstein equations in the spherically symmetric 
model of gravity plus a scalar field 3 '"’. One can use the simplest Schwarzschild coordinates 
(t, r ) which are well-defined for complex r and t, though other coordinate systems may be 
convenient for practical reasons. One starts from wavepackets with small amplitudes 4>o 
which scatter trivially in flat space-time. Then one adds the complex term (2.5), (2.10) 
to the classical action and finds the modified saddle-point solutions. Finally, one increases 
cf >o and obtains saddle-point solutions for the black hole-mediated transitions. The space- 
time manifold, if needed, should be deformed to complex values of coordinates — away 
from the singularities of the solutions. We argued in Sec. 2 that the modified solutions 
are guaranteed to approach flat space-time at t —> +oo and as such, describe scattering. 
The 5-matrix element (1.3) is then related to the saddle-point action S reg in the limit 
of vanishing modification e —> +0. Evaluation of 5-matrix elements is thus reduced to 

33 Another interesting arena for application of the semiclassical method is two-dimensional dilaton grav¬ 
ity [68]. 
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solution of two-dimensional complexified field equations, which can be performed on the 
present-day computers. 

One may be sceptical about the restriction to the spherically symmetric sector which 
leaves out a large portion of the original Hilbert space. In particular, all states contain¬ 
ing gravitons are dropped off because a massless spin-2 particle cannot be in an s-wave. 
Nevertheless, the 5-matrix in this sector is likely to be rich enough to provide valuable 
information about the properties of black hole-mediated scattering. In particular, it is 
sufficient for addressing the questions mentioned in the previous paragraph. 

Furthermore, one can envisage tests of the 5-matrix unitarity purely within the semi- 
classical approach. Indeed, consider the matrix element of the operator 5^5 between two 
coherent states with the mode amplitudes a& and bk (k is the mode wavenumber), 

(a|5^5|6) = f Vc k Vc* k e-f dkc *k Ck {a\S^\c)(c\S\b) , (6.1) 

where on the r.h.s. we inserted the sum over the (over-)complete set of intermediate coherent 
states. For sifficiently distinct semiclassical states | a) and \b) the integrand in Eq. (6.1) is a 
rapidly oscillating function. Then, it is natural to assume that the integral will be saturated 
by a unique saddle-point state |co) which does not coincide with the dominant final states 
of the transitions starting from | a) and |6). This suggests that the amplitudes (a|5^|co) 
and (co|5|6) correspond to rare exclusive processes and can be evaluated semiclassically. 
Substituting them in Eq. (6.1) and comparing the result with the matrix element of the 
unity operator, 

(a| 116) = dka l bk , (6.2) 

one will perform a strong unitarity test for the gravitational 5-matrix. Of course, this 
discussion relies on several speculative assumptions that must be verified. We plan to 
return to this subject in the future. 

Finally, in this paper we have focused on the scattering processes with fixed initial 
and final states which are described in the in-out formalism. In principle, the semiclassical 
approach can be also applied to other quantities, e.g. the results of measurements performed 
by an observer infalling with the collapsing matter. The latter quantitites are naturally 
defined in the in-in formalism with the corresponding modification of the path integral. 
One can evaluate the new path integral using the saddle-point technique. Importantly, the 
new semiclassical solutions need not coincide with those appearing in the calculation of the 
5-matrix. If they turn out to be different, it would imply that the infalling and outside 
observers describe the collapse/evaporation process with different semiclassical geometries. 
This will be an interesting test of the the black hole complementarity principle. We leave 
the study of this exciting topic for future. 
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A A shell of rotating dust particles 

Consider a collection of dust particles uniformly distributed on a sphere. Each partice has 
mass 5m and absolute value 5L of angular momentum. We assume no preferred direction 
in particle velocities, so that their angular momenta sum up to zero. This configuration 
is spherically-symmetric, as well as the collective gravitational field. Since the spherical 
symmetry is preserved in the course of classical evolution, the particles remain distributed 
on the sphere of radius r(r) at any time r forming an infinitely thin shell. 

Each particle is described by the action 



(A.l) 


where in the second equality we substituted the spherically symmetric metric (2.9) and 


introduced the time parameter r. To construct the action for r(r), we integrate out the 


motion of the particle along the angular variable p using conservation of angular momentum 



(A.2) 


5L 


V-gabV a y b - r 2 ip 2 


It would be incorrect to express <fi from this formula and substitute it into Eq. (A.l). To 
preserve the equations of motion, we perform the substitution in the Hamiltonian 


6H = p a y a + 5Lp - 5C 


(A.3) 


where p a and 5L are the canonical momenta for y a and ip, whereas 5C is the Lagrangian 
in Eq. (A.l). Expressing <p from Eq. (A.2), we obtain, 


5H = p a y a + V-Saby a y b VS™ 2 + 5L 2 /r 2 . 


(A.4) 


From this expression one reads off the action for r(r), 



(A.5) 


where we fixed r to be the proper time along the shell. We finally sum up the actions (A.5) 
of individual particles into the shell action 



(A.6) 


where N is the number of particles, m = N5m is their total mass and L = N5L is the 


sum of absolute values of the particles’ angular momenta. We stress that L is not the total 


angular momentum of the shell. The latter is zero because the particles rotate in different 
directions. 
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B Equation of motion for the shell 


In this appendix we derive equation of motion for the model with the action (3.5). We 
start by obtaining expression for the shell energy-momentum tensor. Let us introduce 
coordinates (y a , 9 a ) such that the metric (2.9) is continuous 34 across the shell. Here 9 a , 
a = 2,3 are the spherical angles. Using the identity 

f d 2y ^ 5 {2) (y - y(r)) = 1 , (B.l) 


we recast the shell action (3.2) as an integral over the four-dimensional space—time, 

Ssheii = - J d2 y d ^r J dT m esV-gaby a y b <5 (2) (y - y(r)) • (B.2) 

Here r is regarded as a general time parameter. The energy-momentum tensor of the shell 
is obtained by varying Eq. (B.2) with respect to g a b and r 2 (y), 


rrab _ 2 d S shell _ :a . b m eS f d (2) (y - y{r)) 

T - h ’“ -7^^~ vy ^J ’ 

a = 2 r 2 5S she g = L 2 I" S^(y-y(T)) 
sheila ^ 6r 2 4 7rr 4 m ef f J v/=g 


(B.3a) 

(B.3b) 


where in the final expressions we again set r equal to the proper time. It is straightforward 
to see that the r-integrals in Eqs. (B.3) produce d-functions of the geodesic distance n 
from the shell, 


< 5 (n) 



$ {2) (y ~ y(r)) 


(B.4) 


We finally arrive at 


^shell — tshell ^( n ) 


m e sy a y b 

Airr 2 


+ a _ 
^shell (3 


ra 

b 


L 2 

8tt m e ffr 4 ’ 


(B.5) 


where T shell p oc 5^ due to spherical symmetry. 

Equation of motion for the shell is the consequence of Israel junction conditions which 
follow from the Einstein equations. The latter conditions relate t^ ell to the jump in the 
extrinsic curvature across the shell [56, 57] 


(«!)+ ~ (K)- = - 8 * ( t>‘ - hat 


A 

i v L shell A 


(B.6) 


Here hy is the induced metric on the shell, K^ v is its extrinsic curvature, the subscripts ± 
denote quantities outside (+) and inside (—) the shell. We define both (K /1U )± using the 
outward-pointing normal, n^drX^ 1 > 0. Transforming the metric (3.6) into the continuous 
coordinate system, we obtain, 

(Tfab\ -a -b (K a \ ra a/ 7 ^ iP 71 

)± - -y y /. 9 = > -> (tf- 7 ) 

V r + f± r 

34 Schwarzschild coordinates in Eq. (3.6) are discontinuous at the shell worldsheet. 
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where dot means derivative with respect to r. From Eq. (B.6) we derive the equations, 


r 

r + /;/2 r + f'J 2 _ L 2 | m eS 
yjf 2 + /+ y/r 2 + JZ m eff r 4 r 2 

Only the first equation is independent, since the second is proportional to its time deriva¬ 
tive. We conclude that Einstein equations are fulfilled in the entire space-time provided 
the metrics inside and outside the shell are given by Eqs. (3.6), (3.7) and Eq. (B.8) holds 
at the shell worldsheet. The latter equation is equivalent to Eqs. (3.8), (3.9) from the main 
text. 

The action (3.5) must be also extremized with respect to the shell trajectory y a {r). 
However, the resulting equation is a consequence of Eq. (B.8). Indeed, the shell is described 
by a single coordinate r(r), and its equations of motion are equivalent to conservation of 
the energy-momentum tensor. The latter conservation, however, is ensured by the Einstein 
equations. 


(B.8) 

(B.9) 


C Turning points at L —» 0 


Turning points are zeros of the effective potential V e s, Eq. (3.9). The latter has six zeros 
Vi, i = 1... 6, which can be expressed analytically at L —> 0. We distinguish the cases of 
massive and massless shell, 
a) m / 0: 


r 1,2 = 


m 


2 (M =F m) 


^3,4 =-h 


iL ML 2 \/2ML 5 / 2 
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b) m = 0: 
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e -^/3 
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All turning points approach zero at L —>• 0 except for n 2 in the massive case. Numerically 
tracing their motion as L decreases from L*, we find that the physical turning point A of 
the reflected trajectory is r§ in both cases. 


D Gibbons—Hawking terms at the initial— and final—time hypersurfaces 

Since the space-time is almost flat in the beginning and end of the scattering process, 
one might naively expect that the Gibbons-Hawking terms at t + = ti and i+ = tf are 
vanishingly small. However, this expectation is incorrect. Indeed, it is natural to define 
the initial and final hypersurfaces as 1. 1 _ = const outside of the shell and t- = const inside 
it. Since the metric is discontinuous in the Schwarzschild coordinates, the inner and outer 


- 35 - 




















Figure 13. Final-time hypersurface in the Gaussian normal coordinates associated with the shell. 

parts of the surfaces meet at an angle which gives rise to non-zero extrinsic curvature, see 
Fig. 13. 

For concreteness we focus on the final-time hypersurface. In the Schwarzschild coor¬ 
dinates the normal vectors to its inner and outer parts are 

e = ( 1 /vT-. o, 0, 0), # = (1/Vu, 0,0, 0) . (D.l) 

It is easy to see that the extrinsic curvature K = is zero everywhere except for the 

two-dimensional sphere at the intersection the hypersurface with the shell worldsheet. Let 
us introduce a Gaussian normal frame (r, to, 0 a ) in the vicinity of the shell, see Fig. 13. Here 
r is the proper time on the shell, n is the geodesic distance from it, and 6 a , a = 2,3, are 
the spherical angles. In this frame the metric in the neighborhood of the shell is essentially 
flat; corrections due to nonzero curvature are irrelevant for our discussion. 

To find the components of and p 1 / in Gaussian normal coordinates, we project them 
on r ^ and n^ — tangent and normal vectors of the shell. The latter in the inner and outer 
Schwarzschild coordinates have the form, 

T>1 = ( ^/± ~k ’ ^ °’ °) ’ ^ = + °’ °) ' ( D - 2 ) 

Evaluating the scalar products of (D.l) and (D.2), we find, 

= chV’i — sh. ip± Ttf 1 , shi/}± = —j= ■ (D.3) 

V f± 

As expected, the normals do not coincide at the position of the shell. To compute the 
surface integral in the Gibbons-Hawking term, we regularize the jump by replacing (D.3) 
with 

= ch^(n) — sh^(n) , (D.4) 

where tjj(n) is a smooth function interpolating between ip- and ip+. The expression (3.4) 
takes the form, 

Sgh = j r 2 d 2 0 J dn ^ I< = (ip + - ip -) , (D.5) 

where in the second equality we used ds = dn/ chip for the proper length along the final¬ 
time hypersurface and K = = — chipip' for its extrinsic curvature. Next, we express 
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ip±(r) from the shell equation of motion (3.8) and expand Eq. (D.5) at large r. Keeping 
only non-vanishing terms at r = rj — > +oo, we obtain Eq. (3.24) for the final-time 
Gibbons-Hawking term. 

For the initial-time hypersurface the derivation is the same, the only difference is in 
the sign of £ > ,J ' which is now past-directed. However, this is compensated by the change of 
sign of r. One concludes that the Gibbons-Hawking term at = L is obtained from the 
one at t+ = tf by the substitution r/ —> ri. 

Note that expression (D.5) is valid also in the model of Sec. 4 describing massless shell 
in AdS. It is straightforward to see that in the latter case the Gibbons-Hawking terms 
vanish at r t j- — > oo due to growth of the metric functions (4.3) at large r. 


E Shell self—gravity at order 1/r 


Let us construct the action for a neutral shell in asymptotically flat space-time taking into 
account its self-gravity at order 1/r. To this end we recall that the shell is assembled from 
particles of mass dm, see Appendix A. Every particle moves in the mean field of other 
particles. Thus, a new particle added to the shell changes the action of the system 35 by 

SS = - j Smdr = J dt +(~ <?W 1 ~ v 2 + + y ) , (E.l) 


where v = dr/dt.+ is the shell velocity in the asymptotic coordinates, M is its energy, and 
we expanded the proper time dr up to the first order in 1 /r in the second equality. At the 
leading order in 1/r, 


M = 


m 

\J\ — v 2 ’ 


(E.2) 


where m is the shell mass before adding the particle. Now, we integrate Eq. (E.l) from 
m = 0 to the actual shell mass m and obtain the desired action, 


So 


J dt + ^ — m\J 1 — 


v z + 


m 2 (l + u 2 )\ 
2r(l - v 2 ) / ‘ 


(E.3) 


From this expression one reads off the canonical momentum and energy of the shell, 


mv 2 m 2 v 

y/l - V 2 + r( 1 - v 2 ) 2 


m m 2 (—1 + 4r 2 + v 4 ) 
2r(l — v 2 ) 2 


(E.4) 

(E.5) 


Expressing the shell velocity from Eq. (E.5) and substituting 36 it into Eq. (E.4), we obtain 
Eq. (3.27) from the main text. 

35 Angular motion of the particle gives 1/r 2 contributions to the Lagrangian which are irrelevant in our 
approximation. 

36 In this calculation the 1/r terms are treated as corrections. 
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